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ABSTRACT 

Historical Type la supernovae are a leading candidate for the source of positrons observed through their diffuse 
annihilation emission in the Galaxy. However, search for annihilation emission from individual Type la su- 
pernovae has not been possible before the improved sensitivity of INTEGRAL. The total 511 keV annihilation 
flux from individual SNe la, as well as their contribution to the overall diffuse emission, depends critically on 
the escape fraction of positrons produced in 56 Co decays. Late optical light curves suggest that this fraction 
may be as high as 5%. We searched for positron annihilation radiation from the historical Type la supernova 
SN 1006 using the SPI instrument on INTEGRAL. We did not detect significant 511 keV line emission, with 
a 3er flux upper limit of 0.59 x 10~ 4 photons cm -2 s _I for ~1 Msec exposure time, assuming a FWHM of 2.5 
keV. This upper limit corresponds to a 7.5% escape fraction, 50% higher than the expected 5% escape scenario, 
and rules out the possibility that Type la supernovae produce all of the positrons in the Galaxy (~12% escape 
fraction), if the mean positron lifetime is less than 10 5 years. Future observations with INTEGRAL will provide 
stronger limits on the escape fraction of positrons, the mean positron lifetime, and the contribution of Type la 
supernovae to the overall positron content of the Galaxy. 

Subject headings: ISM:individual (SN1006), supernova remnants, gamma rays:observations 
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1. INTRODUCTION 

Searching for 7-ray lines from supernovae (SNe) remains 
one of the pri mary goals of 7 -ray astrophysics, and also IN- 
TEGRAL ( Winkler et al. 2003), as 7-ray line studies can probe 
the nucleosynthesis and explosion kinematics of these events. 
Many of the nuclear decay chains ( 56 Co, 44 Ti) that produce 
7-ray lines in young SNe also produce positrons, making 
them a prime candidate for the source of positrons for the 
diffuse Galactic annihilation emission. However, there has 
not been a detection of positron annihilation radiation from 
an individual SN or a SN remnant (SNR). The main uncer- 
tainties on the expected annihilation fluxes from young SNe 
are the mean lifetime and escape fraction of positrons. For 
SNe Type la, the lifetime of a positron will be small in the 
initial SN, but may become as high as > 10 4 years as the 
SN expands. The mean positron lifetime will be even longer, 
> 10 5 years, in cases where the thermalization takes place 
in the interstellar medium (ISM), rather t han in the ejecta 
JGuessoum. Ramatv. & Li ngenfelterl ll991l) . If the ejecta's 
magnetic field is weak, and/or radially combed, then 95% of 
the 56 Co decay positrons annihilate promptly , and 5% escape 
to an nihilate on a longer time-scale dChan & LingenfelteJ 
1993). Alternatively, a tangled and strong magnetic field 
would confine ^100% of the 56 Co positrons resulting prompt 
annihilation, leaving a much smaller annihilation radiation 
flux from the delayed 44 Ti decays (r ^85 years). 

Simulations of late optical B & V light curves of Type 
la SNe indicate that the 5% escape of the 56 Co positrons is 
the more probable scenario, ass uming the optical light curve s 
trace the bolometric luminosity (Milne. The & Leising 1999). 
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A later study also showed that 16 of 22 SNe la observed at 
late epochs exhibit the same shape BVRI light curves as the 
initial sampling of 10 S Ne la, stre ngthening the co nclusion 
of 5% escape fraction (Milne. Th e & Leisin g 2001). How- 
ever, excess NIR emission at late epochs was reported re- 
cently from SN 1998bu an d SN 2000cx (Spyromilio et al. 
2004; Sollerman et al. 2004). This excess may be due to emis- 
sion shifting into the NIR wavelength range from the optical 
range. The optical-NIR bolometric light curves for those SNe 
favor full positron trapping. Observations of positron annihi- 
lation radiation from individual SNRs can improve our under- 
standing of magnetic field configurations in these objects, as 
well as determine their role in producing the diffuse Galactic 
annihilation emission. 

The positrons are initially hot, and slow down by 
Coulomb losses, ionization, and excitation. Once they are 
cool enough (less than a few hundred eV), the positrons 
can have charge exchange with neutrals, recombine radia- 
tively with free electrons, and/or annihilate directly with 
free or bound electrons, p roducing a line at 511 keV, 
and a positronium continuum ( Bussard. Ramatv & Drachman 
U979tlGuessoum. Jean & GillardN2005l) . The fraction of the 
photons that are emitted in the line and the continuum depends 
on the properties of the annihilating medium. The recent anal- 
ysis of the SPI data yields positronium fractio ns close to 0.95 
for the diffuse Galact ic annihilation radiation ( Churaz ov et alJ 
120051 Uean et alJl2005l) . For the width of the line, SPI finds 
2.4±0.3 FWH M, if the lin e is approximated as a simple Gaus- 
sian (^hurazov et al. 2005). 6 Note that the positronium frac- 
tion and the FWHM values for the Galaxy are not necessar- 
ily appropriate for individual SNe annihilations. The positro- 
nium fraction depends o n the temperature and the ioniza- 
tion state of the mediu m llGuessoum. Ramatv. & Lingenfelterl 
1 1 99 11 Uean etai]l2005l) . and the FWHM can be broadened by 
Doppler effects if the annihilation takes place in the ejecta. 

6 Recent analysis of the SPI observations of the galactic center region sug- 
gest that the line may be composed of a narrow (FWHM ~ 1.3 keV) and 
broad (FWHM ~ 5.1 keV) component Uean et all2 005). 
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FIG. 1. — SPI significance image in 508.5-513.5 band of the fully coded 
field of view. The ASCA contours of SN 1006 is overlaid to indicate the extent 
and position of the source. The 3° FWHM of the imaging resolution of SPI 
is shown for comparison. The significance scale is plotted on the right. 



For weak and/or radially combed magnetic fields that per- 
mit 5% escape, positron transport simulations estimate that 
roughl y 8 x 10 52 positron s escape from a given SN la 
(Milne /The & Leisindll99 9). When combined with the rate 
of 0.5 Type la SNe per century, a steady state flux of 9 x 10~ 4 
photons cm" 2 s" 1 would be expected if the positrons are as- 
sumed to be annihilated at a distance of 8 kpc (distance to the 
Galactic center). The latest measurements of SPI yields a total 
Galactic 51 1 keV flux of 1.5-2.9 x 10" 3 phot ons cm" 2 s" 1 de- 
pendi ng on the Galactic distribution model ( Knodlsed er et al.l 
2005). The higher end of that flux range is in agreement with 
OSSE on CGR O, SMM, and TG RS observations of 511 keV 
line emission (Milne et al. 2001). The expected flux from the 
5% escape case is therefore 30% - 60% of the total flux, with 
the higher fluxes producing the lower SN fractions. There- 
fore, if positrons escape the ejecta, but are confined to lo- 
cal regions of the SNe, the maps of 5 1 1 should trace the re- 
cent SN la history. The bulge-to-disk flux ratio is confined 
to be within the range of 1-3 (and luminosity ratio of 3-9). 
This ratio is consistent with no Galactic scale diffusion, and 
makes Type la SNe prime candidates for the source of Galac- 
tic bulge positrons, a long with low mass X-ray binaries (see 
Knodlseder et al. 2005 for a thorough discussion of potential 
candidates for diffuse annihilation radiation). 

The most promising individual source candidate to search 
for a positron annihilation line is SN 1006, a recent SN 
thought to be of Type la. Our group observed SN 1006 with 
INTEGRAL for ^1000 ks, with the main goals of character- 
izing the h ard X-ray emi ssion using ISGRI and JEM-X in- 
struments (Kalemci et al. 2005), and searching for 511 keV 
emission line with SPI. In this letter we discuss the results of 
the analysis of the SPI data, and place limits on the 511 keV 
line and positronium continuum emission from SN 1006. 

2. OBSERVATIONS, ANALYSIS AND RESULTS 

SPI is a coded-aperture telescope using an array of 19 
cooled germaniu m detectors for high-resolution spectroscopy 
(Vedr enne et alJl2003l) . It works in 20 keV - 8 MeV band, 
and has an energy resolution of ~2 keV at 5 1 1 keV. The fully 
coded field of view is 16°, and the angular resolution is 3°. 

The INTEGRAL observations of SN 1006 took place in two 
sets. The first 250 ks set was conducted early in the mis- 
sion, between MJD 52650 and MJD 52659 corresponding to 
INTEGRAL revolutions 30 and 32. These observations will 



be denoted as "Set I". The second 750 ks set was conducted 
between MJD 53024 and MJD 53034 during revolutions 155- 
158, which will be denoted as "Set II". For SPI, Set I has all 
19 detectors working, whereas for Set II, the active detectors 
were reduced to 18 after the loss of Detector 2. 

Before any analysis, we filtered out the pointings with high 
Anti-Coincidence Shield rates, mostly occurring during the 
entry and exit of the radiation belts. We used OSA 4.2, 
SPIROS 6. 1, and single-detector events for the analysis. Sev- 
eral background models were tried with "GEDSAT" as the 
main model, which assumes that the background level is pro- 
portional to the product of saturated detector trigger rate and 
live-time. We used the continuum energy band of 523-545 
keV for normalization. We tried normalization with an OFF 
observation (empty field observation in Revolution 130), and 
also using a template file (provided by J. Knodlseder). The 
details of the background methods and i nformation about the 
template file can be found in Knodls eder et al.l (120051) . Finally, 
we tried the mean count modulation method in SPIROS which 
does not require a prior background determination. Maximum 
likelihood method in SPIROS was used for imaging. Note that 
SN 1006 is ^30' in diameter, and is effectively a point source 
for SPI. Set I and Set II were analyzed separately. We used 
5, 7, and 10 keV energy intervals encompassing 511 keV to 
factor in possibilities for both narrow and broad emission. We 
also searched for the positronium continuum in 200-500 keV 
band with SPI. 

Fig-IDshows the SPIROS sigma image in the energy band of 
508.5-513.5 keV, using the mean count modulation method. 
We obtained similar images in larger energy bands, and dif- 
ferent background methods. We have not detected signifi- 
cant 511 keV line emission in any of the energy band inter- 
vals around 511 keV. The positronium component was not 
detected either. Table shows the details of the search and 
3cr sensitivity limits of SPI for the given energy band and ob- 
serving time. The combined upper limit is for 18 detectors in 
Set I and Set II. 



Table 1 
SPI Upper Limits 



511 keVLine 


Set I (250 ks) 


Set II (750 ks) 


Combined 


FWHM (keV) 


Flux" 


Flux 


Flux 


2.5 


1.15 x 10" 4 


0.70 X 10^ 


0.59 X 10^ 


3.5 


1.35 X 10- 4 


0.81 X 10^ 


0.70 X 10^ 


5.0 


1.66 x 10" 4 


0.96 x icr 4 


0.83 X icr 4 


Positronium Continuum 






350-500 keV band 


5.8 x 10" 4 


3.4 x 10^ 


3.0 x 10~ 4 



"All fluxes are 3cr upper limits, ph cm 2 s' 



3. DISCUSSION 

SN 1006 is an ideal candidate to search for 511 keV 
positron annihilation line emission, as it is a well studied, 
historical Type la SN with relatively well determined char- 
acteristics. The distance is ^2.2 kpc, the age is 999 years, the 
angular size is ~30', and the ejecta's current electron density 
is 0.7 cm" 3 llWrnkler. Gupta & LonJ2003HLong et alJl20oi 
Milne !l97lllMiIne. The & Leisingll999HRevnoldsll998l) . 

Predictions of the current 511 keV line and positron- 
ium continuum fluxes from SN 1006 depend upon under- 
standing the complicated interaction of the SN shock with 
the ISM, and the resulting degree of magnetic confine- 
ment of positrons (iRuiz-Lapuente & Spruitl 1 1 9981) . Sim- 
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Mean Lifetime of Positrons fyrs) 

FIG. 2. — Predicted 511 keV line fluxes from SN 1006 as a function of 
positron mean lifetime. The solid and dotted curves are for 5% and 12% 
escape fractions, respectively. The expected SPI (5cr) sensitivity to the 511 
keV line for 3500 ksec, and the current SPI 3<r upper limits for the emission 
from the SNRs are shown for reference. 



ulations have shown that positrons that escape from the 
SN eject a typically have ~ 500 keV of energy upon their 
escap e JChan & Lingenfelterl Il993t iMilne. The & Leisind 
1999). The range of 500 keV positrons in the ISM is such 
that few po sitrons would be expected to thermalize within 
1000 years (Guessoum, Ram atv. & LingenfelterlH99l1) . This 
situation is even worse at the location of SN 1006, where 
the ISM densit y appears even lower th a n fou r media treated 
bv iGu essoum. Ramatv, & Lingenfelter ( 1991). Thus, if the 
shock fails to confine escaping positrons, perhaps even ac- 
celerating the positrons, then SN 1006 would be only a faint 
source of annihilation radiation, even if positrons escaped 
the ejecta. The primary argument against positrons cross- 
ing the shock in this fashion are circumstantial, as the dis- 
tribution of annihilation radiation from SNe would then trace 
the distribution of matter in the Galaxy. It has been shown 
that the diffuse 511 keV em i ssion mapping results from 
OSSE and SPI JPurcell et al] fl997t [Milne. The & Leisind 
120011 iKnodlsed er et all2005t lTeegarde net alJl2005l) indicate 
a strong bulge component, with a bulge-to-disk luminosity 
ratio of 3-9. By contrast, the matter content of the Galaxy 
is conce ntrated in the disk (bulge-to-disk mass ratio is only 
0.3-1.0, R obin et al.l 2003). which appears to eliminate large- 
scale positron transport from source to annihilation site (not 
just from SNe, but in general). 

If the magnetic field in the shock successfully confines the 
positrons, they will repeatedly traverse the inside of the bub- 
ble, being reflected numerous times. Preliminary simulations 
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of this scenario suggest that the positrons will encounter a suf- 
ficient amount of matter to be thermalized and annihilate on 
a 1000 year timescale. However, these simulations are too 
crude for us to confidently predict an annihilation rate at 1000 
years after the SN event. The more detailed simulation will 
be conducted for the analysis of the entire dataset including 
this and the near-future observations. 

In Fig. |2| we show the predicted 511 keV fluxes from SN 
1006 as a function of mean positron lifetime for the 5% escape 
fraction. The figure also shows the current 3a SPI upper limit, 
and the 3a sensitivity level that will be achieved after the ap- 
proved future observations (2.5 Msec). If 5% of the positrons 
indeed escape the ejecta and are locally confined with a life- 
time less than 10 5 years, then SPI is expected to detect 511 
keV annihilation emission from SN 1006 once the additional 
observations are performed. 

Fig. |2] also shows the predicted fluxes for SN 1006 if the 
positron yield is increased by a factor of 2.5, so that SNe 
la could account for all Galactic positrons. The SPI non- 
detection suggests against this possibility if the mean positron 
lifetime is less than 10 5 years. Note that, in case of no escape 
from the remnant, the expected 511 ke V flux from only 44 Ti 
decays is no larger than 4 x 10~ 6 photons cm" 2 s" 1 from SN 
1006, which is far too low to be detectable by SPI. 

We stress that the conclusions above depend strongly on 
the assumption of no escape to the ISM, and thermalization 
in the ejecta. The future observations of SN 1006, and other 
SN la remnants of different ages (both younger and older), 
will provide essential information for understanding positron 
transport in SNRs. A non-detection would still leave an am- 
biguity; the reason can either be ~0% escape fraction from 
the ejecta, or very long thermalization timescale of positrons 
in the ISM. On the other hand, detection of annihilation ra- 
diation would be a strong indicator of thermalization in the 
ejecta, and confirm whether SNe la are important contribu- 
tors of galactic positrons. 
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